Reciprocity in quantum, electromagnetic and other wave scattering 
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The reciprocity principle is that, when an emitted wave gets scattered on an object, the scattering 
transition amplitude does not change if we interchange the source and the detector - in other words, if 
incoming waves are interchanged with appropriate outgoing ones. Reciprocity is sometimes confused 
with time reversal invariance, or with invariance under the rotation that interchanges the location of 
the source and the location of the detector. Actually, reciprocity covers the former as a special case, 
and is fundamentally different from - but can be usefully combined with - the latter. Reciprocity 
can be proved as a theorem in many situations and is found violated in other cases. The paper 
presents a general treatment of reciprocity, discusses important examples, shows applications in 
the field of photon (Mossbauer) scattering, and establishes a fruitful connection with a recently 
developing area of mathematics. 
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1. INTRODUCTION 

Symmetry is a key notion in physics. In the field of 
natural sciences, symmetry can be interpreted as a con- 
cept of balance or patterned self-similarity [lj, regard- 
ing to concrete as well as abstract objects like a phys- 
ical system itself or, for example, a theoretical model, 
respectively @. According to the Oxford online dictio- 
nary, symmetry is a law or operation where a physical 
property or process has an equivalence in two or more 
directions, or events/actions are balanced/equal in some 
way. A synonym to the term symmetry is invariance, 
expressing the fact that a special operation (called sym- 
metry transformation) , which could be a change of some 
physical parameters - like a geometric transformation, a 
change of polarization, parity, charge, the arrow of time, 
etc. - does not change some particular property of the 
system. Depending on the studied system, symmetry 
may have various measures and operational definitions. 
In the theoretical model of quantum mechanics, symme- 
try transformations were given by Wigner Q as general 
operators preserving the modulus of scalar products of 
the vectors of the Hilbcrt space representing the physi- 
cal states. Further, the symmetry theorem of Wigner [3[ 
states that any symmetry transformation can be repre- 
sented by either a unitary linear or an isometric conjugate 
linear (usually called antiunitary) operator. 

A large class of physical processes can be described in 
the framework of scattering theory 0, @ j which is a theo- 
retical sub-model inside quantum mechanics, but can be 
interpreted in classical electrodynamics and in the classi- 
cal mechanics of elastic waves as well. As documented in 
the literature since long ago, physical intuition suggests 
that, when reversing the position of source and detec- 
tor in a wave scattering experiment (see Fig. [T]) , the 
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observed signal will not change. This condition, called 
the reciprocity principle, is indeed fulfilled in most cases. 
Nevertheless, it does not follow directly from first princi- 
ples, therefore, it is not necessarily fulfilled. 

The physical term 'reciprocity' appeared already in the 
19 th century. The first reciprocity equations related re- 
flection and transmission of light at an interface between 
lossless optical media, and were derived by Stokes [f3]. 
The principle of reciprocity was later generalized for more 
complex scattering systems in the field of electromag- 
netic waves 0, Hj], sound waves Q, electric circuits flOf . 
and radio communication [111, as well as in quantum me- 
chanical scattering problems [12| . The reciprocity related 
publications cover the whole 20 th century, as it is summa- 
rized in the review paper of Potton [l3| . Reciprocity can 
be proven for various scattering problems with certain 
limits of validity [IJ-[l7[- 

The number of reciprocity related publications grew 
intensively in the last decade as well, and all the ap- 
proaches that appeared earlier were subject to further 
developments. Among others, the reciprocity relations 
of Stokes @ were generalized to multilayers [IH and to 
absorptive multilayers [l9|, [2(| . The applications in elec- 
trodynamics for electric circuits and antennas were stud- 
ied by Sevgi [2l| , and the reciprocity theorem of classical 
electrodynamics in case of material media containing lin- 
early polarizable and linearly magnetizable substances 
was formulated by Mansuripur and Tsai [22I ]. Nonre- 
ciprocal devices (circulators and isolators) with on-chip 
integration possibility were recently suggested by Kamal 
et al. [23j . In parallel, reciprocity admits applications in 
particle scattering, acoustics, seismology, and the solu- 
tion of inverse problems as well [l3| . 

Reciprocity was considered for nonlocal electrodynam- 
ical (24J . |25| and nonlocal quantum mechanical systems 
[2r| . In a recent publication of Leung and Young [27l |. 
the aspect of gauge invariance was discussed from the 
point of view of quantum mechanical interpretation of 
reciprocity, and new gauge invariant formulations of reci- 
procity were suggested and analyzed. 
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What is reciprocity? In many works it is simply re- 
lated to time reversal symmetry [28|, as was done in 
the well-known reciprocity theorem of Landau and Lif- 
shitz (30j . The optical reciprocity theorem, however, re- 
vealed that absorption, which violates time reversal in- 
variance, conserves reciprocity in polarization indepen- 
dent cases [3l| . which observation was also expressed in 
scattering theory 12j. In parallel, according to the orig- 
inal reciprocity principle, namely, invariance under the 
interchange of source and detector, one could have the 
impression that reciprocity is identical to a rotation by 
180°. However, this latter interpretation also proves false 
since there exist scatterers with no 180° rotational sym- 
metry but fulfilling the reciprocity principle. 

The currently typically used condition of reciprocity 
in linear systems is the self-transpose (also called com- 
plex symmetric) property of the matrix of the scatter- 
ing potential, of the index of refraction, of the dielec- 
tric /magnetic permeability tensors, or of the Green's 
function [3, This condition, however, depends on 
the frame, on the polarization basis chosen. Indeed, ap- 
plying a unitary basis transformation, the self-transpose 
property of the matrix is not conserved, as it can be 
demonstrated on the case of a Hermitian matrix, which 
is not self-transpose in general, but can be diagonalized 
- hence, the self-transpose form is obtained by an ap- 
propriate unitary transformation. In the light of this 
observation, the physical content of reciprocity seems to 
be unclear. Our main task is to give a proper frame- 
independent description of reciprocity, extending the ex- 
cellent early work of Bilhorn et al. [12( • 

The nonreciprocal properties of systems are even less 
understood. Magneto-o ptic al systems are typically cited 
as nonreciprocal media [l3l [23[ , but detailed analyses of 
the reasons of reciprocity violation have not been given. 

Reciprocity violation can be obtained in case of 
magneto-optical gyrotropy which is a well-known 
property of the Mossbauer medium [32]. At Mossbauer 
resonances, the ratio of the time-inversion-violating to 
normal potentials is typically of the order of one thou- 
sand! The resonant Mossbauer medium is absorptive and 
gyrotropic and, accordingly, for well-defined geometrical 
situations, significant reciprocity violation is expected. 

The content of reciprocity is the same for any type of 
classical wave as well as for quantum mechanics. For def- 
initeness, we discuss reciprocity in the quantum mechan- 
ical framework, in scattering theory that corresponds to 
the Schrodinger equation. Note that any classical wave 
equation can be rewritten in the form of a Schrodinger 
equation [as done, for example, by Richtmyer [33J, p347], 
Akhiezer and Berestetskii [34], Ch. I §1J, Taketani and 
Sakata (35{, and Feshbach and Villars [3g]], and, under 
this correspondence, what is probability density in the 
quantum context, is energy density in case of classical 
waves. Especially, the case of two-component wave func- 
tion, which in quantum mechanics describes a 1/2 spin 
particle, e.g., a neutron, is equally able to represent the 
two transversal polarization degrees of freedom of pho- 



ton so the scattering of slow neutrons and that of photons 
admit a common formalism [37], [38]. 

In Part [2J we present the general formalism of 
reciprocity. Part [3] investigates the case of two 
spin/polarization degrees of freedom in detail, and the 
results are illustrated and applied on examples related to 
the area of Mossbauer scattering in Part 2) 

Our discussion analyzes the relationship of reciprocity 
to time reversal invariance and to rotational invariance, 
studies a form of quasireciprocity and the specialties 
emerging in the Born approximation, and uncovers a link 
to a recently expanding area of mathematics the results 
of which assist physics in identifying and finding systems 
with the reciprocity property. 



2. THE GENERAL FORMULATION OF 
RECIPROCITY 

To formulate reciprocity, let us first revisit two top- 
ics involved, scattering theory and antiunitary operators, 
briefly summarizing the ingredients utilized in what fol- 
lows. 



2.1. Notations: Scattering theory 

Concerning scattering theory, we use notations, con- 
ventions, and standard results from Schiff 0] and Mes- 
siah @j_see also Galindo and Pascual |39j and Reed and 
Simon [4fJ, for example, for technical details. 

Let Hq be a self-adjoint Hamiltonian, which, for sim- 
plicity, will be called a free Hamiltonian although it need 
not really describe a free quantum/wave propagation - 
for example, neutrons emitted by a source may travel 
through a guiding magnetic field. Furthermore, let a po- 
tential V describe a scatterer. V is not assumed to be 
self-adjoint, which allows absorption effects to be incor- 
porated. 

With the stationary Green's operators 



:= (E-H±ie)~ 



(1) 



(e \ understood), if u is an eigenstate of Hq with a 
real eigenvalue E then, under suitable conditions on V, 
the states introduced as 



X ± :=u + G±Vu, 
X T± :=u + G%*V*u 



(2) 
(3) 



prove to be such E eigenvalued eigenstates of H = Hq+V 
and its adjoint H\ respectively, that u is their asymp- 
totically incoming ('+' sign) or outgoing ('— ' sign) part. 

For any two eigenstates u a , up of Hq with eigenvalue 
E a = Ep, the u a — > up elastic scattering transition am- 
plitude reads and satisfies 



(f3\T\a) := (u f »V X t) = {Xp-,Vu a ) 



(4) 
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In case the scatterer can be divided into two sub- 
scatterers, V = V\ + V 2 , the transition amplitude can 
also be given as a sum as 

08 \T\ a) = (fi \Ti\a) + {xlp ,V 2X +) , (5) 

where (/3 |Ti| a) is the transition amplitude of scattering 
on V\ alone, and Xib a l so corresponds to V\ only. Nat- 
urally, the role of V\ and V 2 can be interchanged here. 

As for approximations - to which one is forced to resort 
in many applications, - the (1 st ) Born approximation is 
when the scattering solutions are replaced by the corre- 
sponding free ones, 

X + ~ u, x T ~ ~ u, (6) 
and thus ^ is approximated as 

(/3\T\a)tt( U0 ,Vu a ), (7) 

and ([5]) as 

ifi\T\a) « {up,V lUa ) + { Uf3 ,V 2 u a ) . (8) 

We can see that scattering formulae get considerably 
simplified in the Born approximation. 

2.2. Notations: Antiunitary operators 

The key notion behind reciprocity is the notion of 
antiunitary (also called conjugate unitary) operators. 
On a separable complex Hilbert space W, an operator 
U : % — > % is called unitary if it is isometric, 

11^11 = 11^11, ^eH, (9) 

and linear, 

U (A^j + A 2 V 2 ) = A x Z7Vi + A 2 C/V 2; 

A 1 ,A 2 eC, ^6«. (10) 

For antiunitary operators K : % — > % 1 isometry remains 
valid, 

pty|| = Ml, ipeu, (11) 

while linearity is replaced by antilinarity (conjugate lin- 
earity), 

K (A^i + A 2 ^ 2 ) = KKipi + \* 2 Kip 2 , 

A 1; A 2 eC, Ai 2 e«, (12) 

where * denotes complex conjugation. In both cases, 
isometry implies the existence of the inverse operator, 
and the inverse proves to coincide with the adjoint, 

U~ l = U\ K- 1 =K\ (13) 



which are again unitary and antiunitary, respectively. 
With the aid of the so-called polarization identity, from 
the isometric property one finds, for any t/'1j' ! /'2 e H, 

{Urk.,UTh) = {rk. i ik) (14) 

and 

{Kil> 1 ,Kfa) = (ih,il> 1 ), (15) 

respectively. It is this scalar product swapping property 
(Tr5"|) that will be shown below to be the key point why 
reciprocity is connected to antiunitary operators. 

The most frequently treated antiunitary operators are 
the involutive ones, and are called conjugations. Namely, 
an antiunitary operator C is a conjugation if it is involu- 
tive, 

C 2 = I, (16) 

with I denoting the identity operator of the Hilbert space 
H. Similarly defined are the anticonjugations, those an- 
tiunitary operators whose square is —I rather than I (an- 
tiinvolutions) . The most well-known example for a con- 
jugation operator is the standard complex conjugation 

Jip = ip* (17) 

of complex functions ip in an L 2 Hilbert space. Conjuga- 
tions possess various nice properties. For example, any 
conjugation admits an orthonormal eigenbasis ei,e 2 , . . . 
in H with unit eigenvalues, Ci = c 2 = ■ ■ ■ = 1; and, con- 
versely, any orthonormal basis defines a conjugation by 
being its eigenbasis with eigenvalues 1. 

Antiunitary operators are the same in number as uni- 
tary operators, in the standard sense that they can be 
brought into one-to-one correspondence. Indeed, choos- 
ing an arbitrary antiunitary operator - for later purposes, 
let it actually be a conjugation C - any antiunitary K 
can be written in the form 

K = UC, (18) 

where U is unitary. In fact, 

U-KC- 1 (19) 

is a product of two isometric and antilinear operators, 
thus being isometric and linear, i.e., unitary. Conversely, 
the multiplication of any unitary U with C is similarly 
found to give an antiunitary K := UC. 

In quantum mechanics, any symmetry can be given via 
either a unitary or an antiunitary operator. In practice, 
antiunitary cases are much less frequently encountered 
than unitary ones. Two well-known antiunitary symme- 
tries are charge conjugation and time reversal; the former 
being a conjugation (in the above sense, having C 2 = /) 
and the latter being either a conjugation or an anticon- 
jugation, depending on particle number and spin. 

The subsequent considerations involve not only con- 
jugations or anticonjugations but arbitrary antiunitary 
operators. 
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2.3. The reciprocity condition and its consequences 

In a scattering situation as described before, let us as- 
sume that an antiunitary operator A commutes with the 
free Hamiltonian, 

KHoK- 1 = H , (20) 

and also that it connects the potential V with its adjoint 

as 

KVR- 1 = V f . (21) 

Then, for the full Hamiltonian H = Hq + V, we have 

AAA" 1 = A 1 . (22) 

Eqs. ()20j) - (|21[) can be called the reciprocity conditions, 
and A a reciprocity operator for the system - the reason 
for these names will be clear soon. 

A consequence of (|2"D)) is that, if u is an eigenstate of 
Hq with real eigenvalue E, then Ku is also its eigenstate, 
and possesses the same eigenvalue E. 

In parallel, for any real E, (J2"2"|) implies 

A (A - H ±ie)R'- 1 = (E - H^ie) , (23) 
K (E - H ± iey 1 R- 1 = (E-H^Tie)' 1 , (24) 

the latter following from the former. Eq. (j2~4")l formulates 
the reciprocity theorem for the Green's operator (JTJ) : 

KG^R- 1 = Gf. (25) 

Next, we study the consequences of the reciprocity 
property on the scattering quantities. Considering the 
u a — > up elastic scattering transition amplitude, let us 
introduce the notations 

u-a := Ru a , u-fj :— Kup. (26) 

R commutes with Hq and we are treating an elastic pro- 
cess so 

E a =E w =E = E^=:E. (27) 

Applying fl2T} and {25) on the definitions one 
finds 

K X t=X^, K X ±= X p- (28) 

In words, R maps a scattering process to a "reversed" 
one. Hence, for the transition amplitude, we obtain 

(0\T\a) = (u p ,V X +) 

= (RV X +,Kup) = (VtR X +,R Uf} ) 

= (R X i,VRu p ) = (x^-,Vu J ) 

= {a\T\p), (29) 



where ©, (US), flU), dill), ©, and again (gj) have been 
utilized, in turn. 

This result, (|29|) . is the reciprocity theorem for the tran- 
sition amplitude. Why it is in fact the manifestation of 
the reciprocity principle - which has been explained in 
the Introduction - is clear from that it relates a scattering 
process to a "reversed" one. As anticipated in Sect. 12. 2\ 
it is the left-right interchanging property (|15|) that makes 
this "reciprocal" relation possible, which explains why an 
antiunitary transformation is the heart of reciprocity. 

This immediately indicates that reciprocity is not the 
same as rotational invariance under some rotation: Ro- 
tations always act on Hilbert space vectors as unitary, 
not antiunitary, operators. (More on rotations vs. reci- 
procity follows in Sect. 13.51 ) 

Let us observe that the adjoint of (|2"T]) reads 

RV^R' 1 = V, (30) 

which can also be rearranged as 

V 4 = R- X VR. (31) 

Comparing this (|3"Tj) with (|2"Tj) shows that A -1 is also a 
reciprocity operator. In parallel, ()30j) helps us to derive 

R 2 VR~ 2 = R {RVR- 1 ) R- 1 = RV^R- 1 = V. (32) 

This says that V and R 2 commute, or, in other words, 
that R 2 is a symmetry of the system. Naturally, then 
A 4 , A 6 , . . . and A~ 2 , A" -4 , . . . each also commute with 
V. As a consequence, all the antiunitary operators 
A 3 , A 5 , . . . and A -1 , A~ 3 , A -5 , . . . are also reciprocity 
operators of the system. This allows further (121)1) -type 
formulae, in which u a is related not to Ua = Ru a but to 
i%(3) := R 3 u a , etc. 

2.4. Reciprocal partner systems 

If, more generally, A does not fulfill (|21l) but provides 
a connection with another scattering problem with po- 
tential V, 

RVR' 1 =v\ (33) 
then a completely analogous calculation provides 

(J3\T\a)= (xl~ , Vu-^j = (a\T\p) , (34) 

where x an d T correspond to V, according to the sense. 
In such a case we can call the system with H the recip- 
rocal partner of the system with H . 

The relationship between a system and its reciprocal 
partner is a duality type one, i.e., the reciprocal partner 
of a reciprocal partner is the original system. Indeed, 
taking the adjoint of (|3"3")l gives 

AV^A" 1 = V. (35) 



5 



Note that, for any K, the operator defined as 

V := (KVR- 1 )^ = KV^K- 1 (36) 

trivially automatically fulfills (|33p . The nontrivial ques- 
tion here is whether this definition provides just a mere 
abstract Hilbert space operator or a physically reason- 
able scattering potential, for which scattering theory also 
holds. 

In quantum mechanics, if two Hamiltonians H , H are 
connected by a unitary transformation 

H = tyn^ 1 (37) 

with a unitary multiplying operator Uj, = e l< ^'' r ' acting 
on wave functions as 

4>(t,r) = e^'Mt,*), (38) 

then H and H describe the same physical system with 
the wave functions ip mapped to ip according to |38|) 
and any physical quantity operator O mapped to O via 
O = U^OU^ 1 . Eq. is the most general transfor- 
mation freedom under the requirement |?/>| 2 = \ip\ 2 , the 
preservation of the position representation. For a charged 
quantum particle in an electromagnetic field, such a so- 
called gauge transformation of quantities is accompanied 
by the gauge transformation 

V = V - {h/q)dt4>, A = A + {h/q)V<t> (39) 

of the electromagnetic four-potential (q being the charge 
and H the reduced Planck constant). 

If two such gauge equivalent Hamiltonians are con- 
nected by an antiunitary K as 

KHqR- 1 = H Q . KVR- 1 = y f , KHR- 1 = H^ 

(40) 

then this is actually a reciprocity property of one and 
the same physical system, and K can still be called a 
reciprocity operator of this system. It is only that two 
- different but gauge equivalent - representations of this 
system are appearing in the formulae. The gau ge aspect 
has been emphasized by Leung and Young [27| . If there 
is no such gauge connection between a H and a H then 
they do represent two different physical systems that are 
in a reciprocity relationship. 

2.5. Time reversal is a special case 

If a system possesses not only the reciprocity properties 
pO|l -(pT |l but also V = V*,H = H J < then KH = HK so 
K is actually a symmetry of the system. Conversely, if 
H = W and K is a symmetry for both Hq and H, i.e., 
KHq = HqK and KH = HK, then K is a reciprocity 
operator. 

A seminal special case is when K is the time rever- 
sal operator. That is, if there is no absorption and 



time reversal is a symmetry then the reciprocity theo- 
rem (f2"9"| holds and coincides with what is usually called 
the microreversibility property of the scattering ampli- 
tude di. 

Nevertheless, let us observe that, even if time rever- 
sal is not a symmetry - e.g., if absorption is present - 
still we can have a reciprocity theorem offered by the time 
reversal operator, as long as it fulfills the reciprocity con- 
ditions (f2"0]) - (l2"T1) . In such cases, time reversal ensures a 
property which cannot be called microreversibility any 
more but is still a reciprocity property. 

Therefore, reciprocity is not the same as time reversal 
invariance: 

- time reversal is by far not the only possible reciprocity 
operator, and 

- it can be a reciprocity operator (source of a reciprocity 
theorem) irrespective of whether it is a symmetry of the 
system. 

Time reversal is revisited in Sect. 13.41 

2.6. Connection with recent mathematical results 

Similarly to that a Hamiltonian may, but not necessar- 
ily does, admit a symmetry, we may ask how frequently it 
occurs that the reciprocity conditions (|2"0"|) - ([2"T1) are sat- 
isfied by some appropriate antiunitary operator K. Typ- 
ically it is easier to check K against the free Hamiltonian 
and the less easy task is to investigate the validity of (|21l) . 
To get closer to the answer to the latter question, let us 
choose an auxiliary conjugation C arbitrarily, with the 
aid of which we describe the various possible K's by the 
various possible U's via dTBl . 

This way, we can rewrite (|21[) as 

V = (KVK- 1 )' = (UCVC^U- 1 )^ 

= U(CVC- i y TJ- 1 . (41) 

At this point, it is beneficial to recall that, for n x n 
complex matrices - the operators of the Hilbert space C n 

- the adjoint is the transpose of the complex conjugate, 

M t = (M*) T . (42) 

Noting also that, with the aid of the standard complex 
conjugation operator J, which we have already met in 
the context of complex functions [cf. (fl~7|)], we can also 
express M* as 

M* = JMJ- 1 , (43) 
(|4"2"| can be reformulated as 

Af f = (JMJ- 1 ) 7 , (44) 

yielding 

M r =(JMJ~ 1 ) i . (45) 
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This enables us to define the transpose of an operator of 
an arbitrary Hilbert space, with respect to a fixed conju- 
gation operator C, as 



(CVC 



-i\t 



(46) 



This notation makes it possible to re-express the condi- 
tion (gj) as 



v = uv T u-\ 



(47) 



Why this form is worth considering becomes apparent 
when we turn towards the mathematical literature. In- 
deed, there is a recently increasing interest in conjuga- 
tion operators and, more generally, in antiunitary ones 
[4TM47 1 - Especially relevant to our situation is the paper 
by Garcia and Tener [isj . where the authors study ex- 
actly condition (|47[) . in Hilbert spaces C™, i.e., for n x n 
complex matrices V , U . They present a necessary and 
sufficient condition for those V's which fulfill (l47l) with 
some appropriate unitary U ~ in other words, which l/'s 
are unitarily equivalent to their transpose. Here, let us 
only mention some simple special cases. To this end, we 
recall that a complex matrix S is a self-transpose one 
(also called complex symmetric) if 



S T = S. 



(48) 



It is evident that self-transpose matrices are unitarily 
equivalent to their transpose. What turns out is that any 
matrix Z that is unitarily equivalent to a self-transpose 
matrix, 



[UZU~ 



UZif- 1 (for some unitary U), (49) 



also proves to satisfy Z = UZ T U-\ with U := (U T U) ■ 
Similarly can one find that antiskew-self-transpose 
(also called complex antiskewsymmetric) matrices, i.e., 
matrices of the block matrix form 



Bi B 2 
B 3 B4 



B 2 = — B 2 , B 3 = — B 3 , B 4 = B t , 

(50) 

and, more generally, matrices that are unitarily equiva- 
lent to an antiskew-self-transpose one, are also unitarily 
equivalent to their transpose. 

Now, the necessary and sufficient condition of the finite 
dimensional case found in the work by Garcia and Tener 
(45l | can be directly applicable for our infinite dimensional 
Hilbert space situation when, for practical purposes and 
applications, we restrict our search for U (our search for 
K) to some special product form. Actually, it is just such 
a type of restriction that we are going to consider in the 
following sections. In addition, the known finite dimen- 
sional result paves the way for the corresponding infinite 
dimensional theorem to come. In fact, the straightfor- 
ward infinite dimensional extension of the finite dimen- 
sional condition - replacing the building block matrices 
involved by infinite dimensional operators - provides a 



sufficient requirement. It is plausible to expect, but is 
yet to be justified, that this extended condition is not 
only sufficient but necessary as well. 

To summarize, these mathematical results help physics 
to find systems with a reciprocity property. 

As an example, we close this section on the mathe- 
matical literature by quoting a finding by Garcia and 
Tener [45[: For n < 7, any n x n complex matrix that 
is unitarily equivalent to its transpose, V = UV T U~ 1 , 
can be brought into self-transpose form by a certain uni- 
tary transformation, UVU^ 1 = (UVU^ 1 ) 1 . In the light 
of Part [3J it will be apparent that many physical ap- 
plications benefit from this, since the labor of finding a 
reciprocity operator is reduced to finding an orthonor- 
mal basis transformation that makes the potential ma- 
trix/operator self-transpose. 



2.7. Reciprocity violation 

In situations where (P2H)) and (121)1) are valid but ([2"Tj) 
does not hold, it is an interesting question that "to what 
extent" reciprocity is violated. For generic R, this is not 
easy to answer but, for the cases when R 2 commutes 
with V, 



K 2 VR- 2 = V, 



(51) 



it is possible to provide a quantitative solution. For ex- 
ample, K's obeying 



K 2 = e lK I 



(52) 



belong to this class, and ([5"2"j) includes the antiunitary 
operators that occur most frequently in physics, i.e., time 
reversal and charge conjugation. [Actually, this e lK can 
only be ±1, as one finds substituting it into KK 2 = K 2 K 

m 

The advantage to have ([51]) is that it is equivalent to 



(53) 



KVR- 1 = R- l VK 
Therefore, if we decompose V as 

V = V+ + V-, v± 
then ([5"3")l ensures that 



- (V±R- 1 V^K) (54) 



KV±R- 1 =±Vl, 



(55) 



as is easy to check. Taking a look at ([2T]) , we are pleased 
to realize that we have succeeded in decomposing V as 
a sum of a reciprocity preserving term (V+) and a maxi- 
mally reciprocity violating one (V-). Utilizing ([S]) for the 
a — > j5 and j3 — > a amplitudes yields 

(0\T\a) - (P\T+\a) + (xl^V-xt) > ( 56 ) 
(a\T\p) = (a\T + \-p) + ( X ^, K_ X ±) . (57) 
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The component V+ is reciprocity preserving so 

(p\T + \a) = (a\T + \J), (58) 
and thus we arrive at 

(P\T\a) - <S|T|^> = (xl^V-xt) ~ (xl~,V-X±) ■ 

(59) 

Here, it is the right hand side that measures the extent 
to which reciprocity is violated. 

In case we wish to evaluate (|59[) approximately and 
apply the Born approximation based formula (j8]) then 
reciprocity violation is approximated as simply as 

(P\T\a) - (n\T\P) « (u/j.V-tta) - (u m V-Uj). (60) 

After some straightforward algebra on the second term 
of the rhs that is based upon the observation that 

K^VtK = ±vl (61) 

is also implied by ([53")l . (pO)) can be rewritten as 

(/3\T\a)-(a\T\P)^2(u^u a ). (62) 

The effect of the reciprocity violating potential compo- 
nent V- is manifest. 



2.8. Magnitude reciprocity 

In many experiments, one measures magnitudes 
| (f3 \T\a) | rather than the transition amplitudes 
(j3 \T\ a) themselves. Hence, it may occur that reci- 
procity is violated, (/?|T|a) ^ (a\T\(3}, but this is not 
observed because 

\(P\T\a)\ =\(a\T\f)\ (63) 

holds. Situations (|6"3")l can be termed quasireciprocity 
or magnitude reciprocity. Similarly, two systems may be 
magnitude reciprocal partners of each other [a general- 
ization of Sect. 12. 4j . 

A useful connection between (full, true, complete) reci- 
procity and magnitude reciprocity can be established as 
follows. Let us perform a unitary transformation U on 
a system with potential V, by such a U that commutes 
with Hq. Completely analogously to the lines (|2"3"1) - (|2"5|) 
and (|28l) can one obtain 

UGiU- l =Gt, Uxi=xt, UXp=X% (64) 

where is the Green's operator for the system with 

V := UVU" 1 (65) 

and x& e * c - belong to u& := Uu a etc., in system V. Now 
let us assume that U admits u a , up as eigenvectors: 

u& = Uu a = e l5a u a , Ms = Uup = e lSfl up. (66) 



Then we find 

^± = e ««( / + G±y)u a = e i ^Xa ) (67) 
and, proceeding similarly to (|29p . obtain 

(/3\T\a) = (u^VxZ) = (Uu^UVxt) 

= e< S "-^)(l3\f\a). (68) 

We can see that the transition amplitude for the same 
process in the two systems is the same in magnitude. 

Therefore, if not V but such a V is in a reciprocal 
partnership with a V then V is in magnitude reciprocal 
partnership with V. And, reversely, a magnitude reci- 
procity can be converted to reciprocity with an appro- 
priate transformed system. 

This observation will find an important application in 
Sect. 13. 7\ the importance being illustrated in PartlH 



3. RECIPROCITY FOR TWO 
SPIN/POLARIZATION DEGREES OF FREEDOM 

For waves described by a scalar square integrable com- 
plex function, and for multicomponent wave functions 
with a Hamiltonian that is independent of the spin/po- 
larization degree of freedom, the complex conjugation J 
usually satisfies the reciprocity conditions (|20|) and (|2 1 1) . 
The simplest case where the existence of a reciprocity 
theorem is nontrivial is when we have a two-component 
wave function and the Hamiltonian does depend on the 
two-component degree of freedom (in addition to space 
dependence) . 

3.1. Two-component wave functions and the 
reciprocity conditions 

In the rest of the paper, we concentrate on two- 
component wave functions 

*-(£)■ (»> 

where %j>\ , tp2 ar e spatially square integrable complex 
functions. Such wave functions appear in the quantum 
mechanics of a 1/2-spin particle, and the two polariza- 
tions of light can also be described this way (see more on 
this in Part 2]). Let the scattering potential V = V(r) 
be a 2 x 2 matrix valued function so the matrix entries 
Vjk (r) (j, k = 1, 2) are complex valued functions. On the 
other side, let us assume that Hq does not mix the two 
components of the wave function. For simplicity, we take 
the specific choice 

H := -A; (70) 
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in case we describe a 1/2-spin quantum particle of mass 
/x, energy values are considered hereafter rescaled by 
2/j,/h 2 [h 2 /(2n) = 1 convention]. 

A complete set of eigenfunctions of Ho is formed by 
the functions 

u a (r)=p a e^ r (71) 

with eigenvalue 

E a = k 2 a , (72) 

where k Q is an arbitrary wave vector and p a is a vec- 
tor with two complex components p a j (j = 1,2), two 
linearly independent such polarization vectors p being 
considered for any fixed wave vector. Therefore, in this 
concrete setting, the index a comprises the eigenfunction 
identifying quantities as 

a = k a ,p a . (73) 

For describing scattering processes, it is often benefi- 
cial to introduce the retarded (+) and advanced (— ) scat- 
tering amplitudes of the original and the adjoint problem, 
respectively Q: 

/± (k0,a). :=~J e^ r V jk (r) X ± fc (r) d 3 r, 
/^( k/J ,a).:=-i- |e^ k ^(r) X ^(r)d 3 r, (74) 

each scattering amplitude being a two-component quan- 
tity (j — 1,2), and summation over repeated indices un- 
derstood. With these notations, (Q| can be rewritten as 

(P\T\cx)=-4tt (p p J+ (k^a)) 

= -4^(/ T -(-k Q ,/3), Pa ), (75) 

the notation ( , ) used for the scalar product in C 2 as 
well, i.e., (p,q) =p* j q j =V\1i + ??2<?2- 

Consequently, when the reciprocity theorem (|29[) holds 
then it can be expressed via the scattering amplitudes as 

(p p J+(k ,a)) = (k m p)) 

= (/ T -(-k a ,/3),p Q ) = (f T -(~k^a),pp). (76) 

the latter line following from (|75|) . 

As has been mentioned in Sect. 12.21 the antiunitary 
operators are very many in number. In what follows, 
we will restrict our attention to an important class of 
antiunitary operators, namely, to those K which are of 
the form 

K = UJ, (77) 

where U is a 2 x 2 unitary matrix - mixing the two com- 
ponents of our two-component wave functions, in a space- 
independent way - and J is the antiunitary operator of 



complex conjugation [cf. (flTl) ]. So to say, these are essen- 
tially 2x2 antiunitary operators. With such a K, (j2"B"|) 
gets concretized as 

u w (r) = (Ku a )(r) = (Up* a )e-^ r , 

ka = -k Q , p-=Up* a , (78) 

which can also be written as 

a=-k a ,Up* a (79) 

[cf. m. 

Clearly, a K of this specialized form commutes with 
Hq = —A so the first of the reciprocity requirements, 
(f2"U|). is fulfilled. In parallel, the other - and much more 
nontrivial - reciprocity condition, (|21[) , simplifies to 

V = UV T U-\ (80) 

as we have already seen at (l47l) . Note that, from the 
point of view of Sect. 12. 6[ now C = J and V T is directly 
the standard matrix transpose of V. 

Before starting to analyze condition , we make two 
remarks in passing. The first is that, should V be a non- 
local potential, i.e., an operator acting on wave functions 
as 

(V^)(r) = J ^(r,rXr')dV, ( 81 ) 

requirement (|80[) is concretized as 

V(r,r') = UV(r',r)U-\ (82) 

since 

VHr,r') = V(r',r)*. (83) 

Reciprocity for nonlocal potentials has been considered 
in HIH. 

The second remark is that a larger family of possible 
reciprocity operators is also allowed: Those when U of 
(1771) is a Hilbert space operator acting on wave functions 
as 

(Ety)(r) = W (O^r) (84) 

with a unitary 2x2 matrix Q and an orthogonal trans- 
formation O. Such more general Ks also commute with 
Hq = —A, and allow for a broader range of applications. 
The thorough discussion of such operators is, however, 
more complicated so here we focus on the case of purely 
2x2 matrix Us. 



3.2. Unitary equivalence of a two-by-two matrix to 
its transpose 

Condition (|80|) says that 

V(r) = UV T {r)U- 1 (85) 
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has to be fulfilled at any location r, with a certain 2 x 
2 unitary matrix U. As a first step in analyzing this 
requirement, let us consider it at a fixed r. fn other 
words, let us first pretend that our potential is a constant. 

Then the question is that which 2x2 matrices V are 
unitarily equivalent to their transpose V T . Let us now 
answer this question. 

A useful characterization of any 2x2 matrix V is done 
by four complex numbers Vo,Vx,V2,v 3 via the following 
expansion: 



V = v a Q + v 1 a 1 + v 2 <J2 + v 3 a 3 



VqCTq 



Vq + V 3 Vi - 
Vi + iv 2 Vq 



IV 2 
- V 3 



(86) 



with Vq = ^trV and v = itr(a"V), where the matrices 



<T 2 



1 
1 

-i 

1 



CT3 



1 

1 

1 

-1 



(87) 



i.e., the 2x2 identity matrix and the Pauli matrices, 
are each self-adjoint. Here, v is frequently referred to 
as the Poincare vector corresponding to V. In the case 
of a self-adjoint V, all the coefficients vo, . . . v 3 are real, 
and, thence, v is a real three-dimensional vector, a vec- 
tor in R 3 . When V describes absorption, too, then the 



since ao, arL d a 3 are self-transpose and a 2 is anti-self- 
transpose. Geometrically, transposition means a reflec- 
tion of both Re v and Im v in R 3 , with respect to the 
<J\-o 3 plane. In notation, V i— > V T means 

Rev i-^ Pi 3 (Rev) , Im v h-> P i3 (Imv) , (93) 

P13 standing for the reflection in question. 

To summarize, the problem whether a V admits a U 
with which V — UV T U^ 1 is translated to the geometric 
question whether a rotation exists that brings P13 (Re v) 
back to Rev and P 13 (Imv) back to Imv. 

For the answer, first let us observe that, as has been 
mentioned above, a rotation rotates the component or- 
thogonal to the axis, in the plane orthogonal to the axis 
and preserves the component parallel to the axis, thus the 
difference O n x is orthogonal to n (lies in the plane 
orthogonal to n). Since the difference Rev — P13 (Rev) 
is perpendicular to the a±-a 3 plane, the axis of any rota- 
tion that brings P13 (Rev) to Rev must be within the 
CT1-03 plane. 

Now, a rotation O n (p rotates any plane that contains 
n to another plane containing n, the included angle be- 
tween the two planes being (p. 

Consequently, to any n within the o\-a 3 plane, let us 
consider the plane spanned by n and P13 (Rev) (which 
is unique as long as Rev 7^ 0). Let 1? denote the an- 
gle between this plane and the a\— 03 plane. Then the 



imaginary part of v, again a vector in R 3 , is nonzero. rotation O n ip with tp := 2z9 brings the plane of n and 



Furthermore, for a unitary U, this decomposition can be 
shown - e.g., using [f| p546] - to read 



U 



e lS (cos |cr 



1 sin -ncr 

2 



< S, ip < 2tt, 



where n is a real unit three- vector (element of R 3 ). 

Now, the key observation for our purposes is that, if 
X is a matrix with real components (x Q ,xi,x 2 ,x 3 ) = 
(x ,x), then 

UXU^ 1 = U (x a + xer) U~ l = x a + (O n v x) cr, 

(89) 

where O n v is the rotation of real three-component vec- 
tors around n by angle ip. both defined by ([88]) . Eq. (|89|) 
can be proven applying [5( 

(acr)(bcr) = (ab)cr + i(a x h)cr (90) 

and the geometrically easy-to-check formula 

O n i(1 x = (nx) n+ [x — (nx) n] cos ip + (n x x) sin^ (91) 

[telling that a rotation does not change the component 
parallel to the axis (first term), and rotates the compo- 
nent orthogonal to the axis, in the plane orthogonal to 
the axis (second and third terms)]. 

Being ready to turn towards the question of V — 
UV T U~ 1 , let us observe that 



V T = (v a ao + Viai + v 2 a 2 + v 3 a 3 ) 
= voao + viai - v 2 a 2 + v 3 a 3 , 



(92) 



P13 (Rev) to its <7i-<73-reflected, i.e., the plane of n and 
Rev , and brings no other plane to its CTi-173-reflected. 

Hence, if we want P13 (Imv) to be rotated to Imv, 
too, then P13 (Imv) must be contained in the plane of 
n and P13 (Re v) . We can always choose such an n that 
this is satisfied: 

(1) If the plane spanned by P13 (Rev) and P13 (Imv) 
intersects with the o\-u 3 plane at a line then let n be 
along this line; 

(2) if Rev and Imv arc the multiples of each other, out- 
side the a\-a 3 plane, then any n within the o\-o 3 plane 
suffices; 

(3) at last, if both Rev and Imv are within the <J\-u 3 
plane then no rotation is needed (we choose the identity 
transformation) . 

We have concluded that any V admits a U that can 
be chosen based on Rev and Imv, the rotation O ntp 
corresponding to U being determined uniquely if Rev 
and Imv are linearly independent. 

Having answered our problem, we can make the fol- 
lowing useful 'by-product' observation. If we rotate by 
— d instead of 2d then we rotate both Rev and Imv into 
the <7i— cr 3 plane. In other words, we establish a unitary 
transformation of V into a self-transpose matrix. 

The fact that any 2x2 matrix is unitarily equivalent 
to a self-transpose one has been known for a long time 
- see the mathematical literature invoked in Sect. 12.61 - 
and here it emerges in a geometric incarnation. 

Furthermore, it is adequate to recall again the results 
mentioned in Sect. 12.61 revealing that, in dimensions 
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n < 7, unitary equivalence to the transpose holds for the 
same matrices as unitary equivalence to a self-transpose 
matrix. 

During the above consideration, the special case when 
Rev and Imv are the multiples of each other behaved 
in a distinguished way. This occurs when, in the de- 
composition (|86p. v = cb with a complex c and a real 
three-vector b: 

V = v a + cber. (94) 

These V can be called univectorial, and the other ones, 
with linearly independent Rev and Imv, bivectorial. 
The family of univectorial potentials is important from 
the physical point of view, as it includes the scattering 
interaction of a neutron with an absorptive medium with 
a magnetic field - where b is the magnetic field B it- 
self - and also the Mossbauer scattering of photons (see 
Part [4J - where b is not directly the magnetic field but 
is determined by it. Univectorial potentials will play two 
important roles in our subsequent discussion, one uncov- 
ered in Sect. 13.71 and the other, based on their special 
property in the above geometric picture, turning out in 
the following section. 

3.3. A criterion ensuring a reciprocity operator 

Now, we are prepared to make the second step, taking 
the space dependence of V into consideration. The cor- 
responding second question is: When does (|85[) hold for 
any location r? 

Fortunately, the consideration we made in the above 
section enables us to give the answer immediately. 
Namely, any given U means a given plane - which it 
can rotate to its reflected - so if and only if all Rev(r), 
Imv(r) of all locations r are in a common plane within 
M 3 then a U we seek exists, and it is actually the one cor- 
responding to this common plane. Potentials with such 
a common plane may be called globally uniplanar. 

One remarkable special case is a potential that is glob- 
ally (everywhere) univectorial [see (|94|) ] . For such po- 
tentials, all those locally single vectors b(r) have to fall 
within a common plane. 

Another example, important for applications like those 
considered in Part [4J is when the potential is piecewise 
constant, taking a value V\ on one space domain, a value 
V% on another domain, etc. Then the vectors 

Revi, Imvi, Rev2, Imv2, Rev„, Imv„ (95) 

must be within a common plane. Naturally, here again, 
rotating by half the angle (see the previous section) 
means a unitary transformation of each of these poten- 
tials to self-transpose ones. 

For univectorial potential values V\ , Va , ■ ■ ■ , V n with 
bj, b 2 , . . . , b„, this criterion of a common plane implies 
that any two univectorial values admit a common uni- 
tary transformation to their transpose (and another one 



to self-transpose ones), the common plane being spanned 
by hx and b 2 , being orthogonal to hi x b 2 , with the ex- 
ception of collinear bi, b 2 when the common plane is not 
even unique. The existence of a common plane becomes 
nontrivial only for n > 3. 

We close this section by mentioning the relationship 
between commutativity and joint reciprocity of 2 x 2 ma- 
trices. This question is a natural one in the light of the 
well-known connection between commutativity and si- 
multaneous diagonalizability of diagonalizable matrices. 
Now, omitting straightforward details, one can find that 
V1V2 = V2V1 if and only if 

Re vi x Re v 2 — Im vi x Im v 2 

= Revi x Imv 2 + Imvi x Rev 2 = 0. (96) 

It follows that Re vi , Re v 2 , Im Vi and Im v 2 all lie within 
a plane. Hence, V\ and V2 admit a joint reciprocity op- 
erator. The converse direction does not hold: For exam- 
ple, V\ — o\ and V 2 = 03 do not commute but share a 
common U that connects both of them with their own 
transpose. In fact, U = (J?) suffices as a\ and 03 are 
already self-transpose. Consequently, commutativity im- 
plies, but does not follow from, joint reciprocity. 

3.4. A special case of reciprocity: Time reversal 

Similarly to how KVK^ 1 = V' is simplified, via the 
decomposition K = UJ, to V — t/U 7 ?/ -1 , the gen- 
eralization KVK^ 1 = seen at ff33|) becomes V — 
UV T U~ 1 . A physically important application of this re- 
ciprocal partnership between different systems, discussed 
in Sect. 12. 3[ occurs for Ut = — i&2, when K T — UtJ is 
the time reversal operator of 1/2-spin quantum mechan- 

icsSg. 

Using (|90|) and (|92|) . it is easy to check that the time 
reversal Kt maps 

V = noco + V\(Ji + V202 + (97) 

to 

V = U T V T U^ 1 — 02 (vqctq + — V202 + V303) er^ 1 
= v a - v\cx\ - i; 2 er 2 - 
= vqUq — ver. (98) 

Two remarkable specializations are 

(1) when v = 0, and 

(2) when v = 0. 

In case (1), we find V = V, which means that, for 
spin/polarization-independent potentials, time reversal is 
a reciprocity operator. It is essential to emphasize that vq 
is allowed be any complex number for this so reciprocity 
holds true also in cases when time reversal is not a sym- 
metry (when Imuo 7^ 0). This situation serves as an 
example for what Sect. [2~51 has pointed out: Reciprocity 
is not the same as time reversal invariance. 
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In case (2), we have V = — V, that is, time reversal has 
mapped the potential to its negative. When V — cBcr - 
magnetic dipole in a magnetic field - this means B = B. 
This is in conformity with the well-known fact that time 
reversal must be accompanied by the reversal of mag- 
netic field to obtain a generalized symmetry (a spectrum 
preserving equivalence between systems) , which is now a 
reciprocal partnership as well. Reciprocity is more gen- 
eral than time reversal invariance. 

For clarity, let us remark that the change B = B 
cannot be covered by a gauge transformation (|3"9")) . Gauge 
transformations of the electromagnetic potentials always 
keep the electric and magnetic fields invariant. B = — B 
does indeed mean two different physical systems. 

3.5. Reciprocity vs. rotation 

The physical setting for testing the reciprocity theorem 
(P\T\a) = (a\T\~/3) is the following: A source emits a 
wave from direction k Q with polarization p a towards the 
scattering object, and the scattered wave is detected in 
direction kp, with polarization pp; and this scattering 
is compared to when the incoming direction is kg- with 
incoming polarization p-p-, and the detector direction is 
k^-, detecting polarization p^- Now, in experiments at 
large facilities it would be difficult to change the direction 
of the source. It is much more feasible to perform a 
rotation of the small scatterer sample. Let us now derive 
a version of the reciprocity theorem that does not require 
to modify the incoming direction, which can be achieved 
with the help of an appropriate rotation. 

Let us start by repeating that, by (f79")). we have now 

k« = -k Q , kp = -kp. (99) 

Recall also that we are considering elastic scattering pro- 
cesses (see Sect. l2.1| ). which ensures E a = Ep and implies 
k a = kp [cf. (|72p ]. This means that the vectors k Q and 
kg can be interchanged by some rotation. More closely, 
we are interested in a rotation that brings kg = —kg 
to k a . Actually, one can observe that the rotation by 
7r = 180° around the direction of the momentum trans- 
fer k,g — k Q realizes this desire and, at the same time, 
maps k„ = —k a to k^, which has the benefit that the de- 
tector direction also remains the same as for the process 
a 

Furthermore, in the special case of forward scattering, 
k^ = k Q , when there is no momentum transfer, any di- 
rection orthogonal to k a suffices for the same purpose. 

Let Or denote this rotation around this direction and 
angle, 

n R = (kp - k a ) / |kg - k a \ , P R = n (100) 

(except for kp = k^ , when is arbitrary up to ur k a = 
0). It is known from 1/2-spin quantum mechanics that 
a rotation is represented as a Hilbert space operator on 
the wave functions as 

(ity) (r) = U R i, (O^r) , (101) 



where Ur is the SU(2) transformation [see also ((88)) ] 

TT P R ■ ■ PR ■ fir>r>\ 

U R = cos—a -tsm—n R a=-in R a. (102) 
Especially, 

(Jfcfc) (r) = U R p^°^ = U R Up* a e"^ r (103) 

[note k-(0^ 1 r) = (0 R k) -r], giving for the eigenfunction 
index [cf. and ([79])] 

Ra = R k^ U rPw = -kp, U R Up* a . (104) 

Our H is rotation invariant so, with V R = RVR -1 , 

RHR- 1 = H + RVR- 1 = H Q + V R , (105) 
RGtR- 1 = (G±) R) (106) 

where (G^) R denotes the Green's operator corresponding 
to V R - see (H|) with Hamiltonian Hq + Vr. This enables 
us to derive, from ([2]), 

R X ± = Ru^+RG^R- 1 (RVR- 1 ) Ru- 

= u m +{Gt) R u Rs . (107) 

Consequently, the reciprocity theorem can be re- 
expressed first as 

(P\T\a) = (a\T\p) = (ug, V X ±) = (Rus, RVBT^Rx*) , 

(108) 

and then in final form, 

(kp,pp\T\k a , Pa ) = (kp,URUp* a \T R \k a ,URUp*p), (109) 

where the eigenfunction indices have been explicitly dis- 
played, and Tr stands for the transition amplitude for 
scattering on Vr. 

One can observe that the rotation is able to transform 
the incoming momentum to the outgoing one and vice 
versa, but unable to transform an incoming polarization 
to an outgoing one (and vice versa). Reciprocity is not 
the same as a rotational invariance. 

It can be practical to choose the z coordinate axis of 
our coordinate system parallel to kg — k Q . Then = 
(0 1) T , Ur = —i(J3- Since the transition amplitude, 
being a Hilbert space scalar product, remains invariant 
if we multiply both the initial and the final state by i, 
([TU9"|> is simplified to 

(kp,pp\T\k a ,p a ) = (kp,a 3 Up* a \T R \k a ,a 3 Up*p). (110) 

In the special case U = (o ?), displaying explicitly the 
two components of the polarization vectors, 

(kp,PPi,PP2 \T\k a ,p al ,p a2 ) 

= (kg,p a *,-p Q 2 \T R \k a ,ppl,-ppl) ■ ( U1 ) 
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On the other side, in forward scattering, = k a , a 
distinguished choice for the z axis is the direction of k Q . 
If iir is in the x direction then Ur — —i(T\, and the 
analogous result is 

(k ,p \T\k a ,p a ) = (k^aiUpllTR^cnUp}), (112) 

which, for U = ( J 5 ) , simplifies to 

(kp,pp v pp 2 \T\ k a ,p al ,p a2 ) 

= (kp,Pa2,Pa* \ T r\ k-cPp^Ppl) . (113) 

Similarly, if is in the y direction then Ur = —icr 2 , 
leading to 

{k/3,p/3\T\k a ,p a ) = (kp,-ia 2 Up a \T R \k a ,-ia 2 Up*p), 

(114) 

and, for C/= (I?), to 

(k/3,P/3 1 ,P/3 2 \ T \k a ,p al ,p a2 ) 

= (k/j, -p a2 ,p a $ \T R \ k a , -pp* 2 ,ppl) . (115) 

For reciprocal partners V and V = UV T U~ 1 , the gen- 
eralization of (|109|) reads, naturally, 

(kp,pp\T\k a ,p a ) = {kp,URUp* a \T R \k a ,URUp*p), (116) 

where T corresponds to V and Tr to the rotated RVR -1 . 

If the scattering potential V is spin/polarization inde- 
pendent, V(r) = uo(r)(7o, then not only the free solu- 
tions u a (r) are of the form (|7T|) but the corresponding 
scattering solutions Xa( r ): Xa ± ( r ) are al so a product of 
the polarization term p a and a space dependent func- 
tion, the latter being the appropriate scattering solution 
of the reduced problem of scalar wave scattering on po- 
tential vo(r). Accordingly, the transition amplitude is 
also factorizable, as 

(k , Pf} \T\k a ,p a ) = {p ,p a ) (k^|T red |k a >. (117) 

As a consequence, the reciprocity theorem combined with 
rotation, viz., (| 109[) is easily seen to get simplified to 

(pp, Pa ) (kp\T™ d \k a ) = (p ,p a ) (k p \T R cd \k a ) (118) 

[as (U R Up* a ,U R Up*p) = {p* a ,p*p) = {pp,Pa)}, where the 
rhs of (|118l) refers to the reduced scalar scattering on the 
rotated scalar potential 

M B (r) = ^(O^r). (119) 

This shows that, for polarization independent scattering 
- but only in those cases - reciprocity acts the same way 
as a rotation. 

Historically, reciprocity was first studied for polariza- 
tion independent phenomena. This explains why it had 
to be a later step to recognize that, in polarization de- 
pendent wave scattering, reciprocity deviates from a ro- 
tation invariance, the difference manifesting itself in the 



polarization degree of freedom. Naturally, the third step, 
formulating reciprocity for general quantum/wave sys- 
tems, as done here in Part [3J makes it apparent that 
reciprocity is related to an antiunitary operator and ro- 
tations to unitary ones. This remarkable mathematical 
difference carries considerably different physical content. 

3.6. Fixing the reciprocity operator to processes 

Reciprocity and its violation has already been dis- 
cussed at the general level, in Part [2] It is interesting 
to observe that, for two-component wave functions, there 
is a special additional possibility. Namely, if we have a 
Hamiltonian Hq + V and are interested in the relationship 
between any two scattering processes a — > j3 and /3 — > a 
then there may be a if, unique in the range chosen at 
([77]). which satisfies 

u-a=Ku a , u-p=Kup. (120) 

In fact, using the decomposition (J77J, the two formulae 

Pa = u P * a , Pp = u P ; (121) 

define a linear operator U uniquely - as long as p a and pp 
are linearly independent - and if this U is unitary then 
we have arrived at a uniquely defined antiunitary K. 

Then we are allowed to ask whether V = UV T U~ 1 , 
and if not then how much reciprocity is violated, along 
the lines of the general treatment of Sect. 12.71 

This is a reverse approach in the sense that not u a , up, 
and K define «„ and u-g, as we proceeded at (|26j) . but u a , 
up, Ua, and u-e define K. This special possibility, which 
exists only for two polarization degrees of freedom, en- 
larges the range of application of the reciprocity theorem 
compared to the general case. 

Actually, a very frequent approach is to fix processes 
- momenta and polarizations - and to ask whether the 
corresponding K is a reciprocity operator for the system's 
V. In fact, typically (tacitly) U = <tq, K = J is assumed 
and thus the question is the self-transposeness of V. 

Now, for n < 7 polarization degrees of freedom, 
indeed the reciprocity property is equivalent to self- 
transposeness with respect to an appropriate orthogo- 
nal basis (cf. Sect. 12. 6p . However, it is beneficial to 
be prepared for other polarization bases. Furthermore, 
there are physical situations with dimensions n > 7, like 
many-particle quantum mechanics of spin-half particles, 
or the presently considered two-component wave situa- 
tions when we do not restrict ourselves to 2 x 2 matrices U 
but allow U be some essentially infinite dimensional uni- 
tary operator of the infinite dimensional Hilbert space. 

That being said, the approach to ask 2x2 self- 
transposeness fits numerous physical settings and is a 
convenient and simple means of finding situations pos- 
sessing reciprocity. 

As a matter of fact, an extended version is to choose 
two orthogonal polarizations p a ,Pb and to consider four, 
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rather than two, processes to relate: 

Pa -> Pa, Pa ~> Pb, Pb ~> Pa, Pb ~> Pb- (122) 

(Some k Q ,k^ must also be fixed, and a rotation of the 
previous Section may also be included.) If 



/<-. = !„]. I? 



(123) 



then the question of reciprocity is the question of whether 
Vjk = Vkj (j, k = 1,2). Among these, the only nontrivial 
one is whether or not 



V 12 = V 21 



(124) 



holds. This single condition ensures reciprocity for four 
processes. 



3.7. Potentials exhibiting magnitude reciprocity 

As has been raised in Sect. 12.81 many experiments mea- 
sure only the magnitude of the transition amplitude, and 
one can speak of magnitude reciprocity when two tran- 
sition amplitudes, related by an antiunitary K, are the 
same in magnitude. There, we derived a possibility for 
this to occur, which consideration, actualized for the sit- 
uation (|122[) - (1124p . implies that, if V is self-transpose 
and 



u = 





a i5 b 



then 



v ■= uvu- 1 



(125) 



(126) 



exhibits magnitude reciprocity for K = J and processes 
In detail, (|126[) reads 



Vii V 12 
V 2 i V22 



Vu 



so V12 = V 2 \ implies 



V 12 = e ld V 21 



(127) 



(128) 



with some angle S. This means that, similarly to that a 
self-transpose potential is reciprocal with respect to K — 
J, a "phase self-transpose potential" [i.e., a potential 
obeying (I128p ] is magnitude reciprocal with respect to 
K = J. 

An example for a phase self-transpose potential is an 
everywhere univectorial potential [cf. (|94pj. i.e., 



with real b(r) and possibly complex vo(r),c(r), if it obeys 
some restriction. That restriction can be revealed via the 
explicit matrix form of such a potential, 



V = 



Vq - 



cb 3 c (61 

h ib 2 ) Vq - 



-ib 2 ) 
cb 3 



Vq -I- CU 3 

cb l2 e lSl2 



cb 3 cb 12 e 

V - Cb 3 

(130) 

with writing b\ +ib 2 in the polar form b\ 2 e lSl2 . From this 
we can read off that, if 5i 2 is space-independent, then this 
potential is phase self-transpose. 

We close this Section with the simple observation - 
which is easy to derive from ([89)1 - that a conjugation 
(|126p by a U (|125|) acts on V like the unitary transfor- 
mation of the polarization basis (|123j) that corresponds 
to some rotation around the a 3 axis. This side remark 
will find application in Part |4j 



3.8. Forward transmission processes 

There is an important special type of forward scatter- 
ing (k^ = k a = k, k = |k|) settings where the scattering 
solution \a [cf. ©] equals u a before reaching the scat- 
terer object and is TfU a afterwards, where the 2x2 for- 
ward transmission matrix Tf is an analytic function of 
V. For example, in Mossbauer optics, the wave crosses a 
layer of width d perpendicularly, the potential V is con- 
stant within the layer and is zero outside, and 



Tr 



j {kd+ £ V ) =goao + gV} 



(131) 



where the exponential has been expanded using the iden- 
tity 



sin Vdet N 



e M = e |trM cos Vdet Na + — N 



VdetN 



N = M — - (tiM) cr 



(132) 



for 2x2 matrices M - a straightforward consequence 
of ((90)) - and the V dependence of the complex multi- 
plier scalars g$, g is not displayed because only the matrix 
structure of Tf will be relevant for the present consider- 
ations. 

Indeed, from this expanded form it is easy to observe 
that T f is 

- self-transpose, 

- phase self-transpose, 

- univectorial 

if and only if V is respectively 

- self-transpose, 

- phase self-transpose, 

- univectorial. 

Now, for such forward transmission situations, the task 
to calculate (/? \T\ a) is reduced to determine (pp,Tfp a ). 
In accordance with this, the reciprocity theorem of the 
transition amplitude simplifies to 



V(r) = v (r)a Q + c(r)b(r)er 



(129) 



(Pt3,T f p a ) = (pa,T f Pp). 



(133) 



14 



Further, if V is univectorial then 612 [cf. (I130[) ] is space- 
independent so V is phase self-transpose. However, more 
is true for such cases. Namely, a unitary transformation 
of the polarization basis (|123l) acts on b of a univectorial 
V, hence, S12 of a space-independent b remains space- 
independent after any rotation of b. The consequence is 
that, in these forward transmission settings, a homoge- 
neous univectorial V exhibits magnitude reciprocity with 
respect to any (essentially 2x2) antiunitary operator K. 

The question of reciprocity for transmission through 
more than one such layer can be analyzed the same way 
as for general forward scattering in the Born approxima- 
tion (see the following Section). 



3.9. The Born approximation 

The importance of the (1 st ) Born approximation [cf. 
©~©] nes in its validity for many situations, e.g., for 
scattering on thin enough layers. In parallel, from the 
aspect of reciprocity [cf. (f60|) and ([62]) ]. it plays a spe- 
cial role as the violation of reciprocity or magnitude reci- 
procity may vanish in the Born approximation. Let us 
now consider its details for two polarization degrees of 
freedom. 

In the Born approximation, the scattering amplitudes 
(|74[) become linear in the incoming polarization: 



f(kp,a) j « f(kp,k a ) jk p k 



(134) 



for each of the four types of scattering amplitude. Cor- 
respondingly, the transition amplitude simplifies to the 
two-component scalar product 



(pp, [J e l ( k °- k ^ r F(r)d 3 r 



Pa 



In case of forward scattering, it reduces to 



Pa 



(135) 



(136) 



If V is piecewise constant, Vi on a spatial region of volume 
V/, with I = 1, ...,n {e.g., a sample constituted by n 
homogeneous layers) then (|136p gives 



(137) 



If, further, each Vi is univectorial such that v; = cb/ with 
/-independent c then the sum is a univectorial matrix 
with Poincare three- vector V;b/ and, consequently, 
we have magnitude reciprocity in any orthonormal polar- 
ization basis [see (13.71)1. 



4. EXAMPLES AND APPLICATIONS 

The general treatment of reciprocity (and nonreciproc- 
ity) given in Part [5] has been investigated in large detail 
for two-component wave functions in Part [3J because of 



the many related important applications in elastic scat- 
tering of photons and neutrons. 

We recall the result given in Sect. 13.91 that nonrecipro- 
cal forward scattering disappears under the conditions of 
the first Born approximation, namely, in case of weak 
scattering. The recoil-less nuclear resonance forward 
scattering of photons, known as Mossbauer scattering, 
can realize the case of strong scattering. Indeed, soon 
after the discovery of the Mossbauer effect the so-called 
blackness effects were reported and identified as the re- 
sult of multiple scattering. In the theoretical description 
of the Mossbauer scattering given by Blume and Kistner 
[32| . the ensemble of nuclei and electrons, as scattering 
centers, represent an anisotropic and absorbent optical 
medium, which is described by a 2x2 complex index of re- 
fraction n :37j corresponding to the two possible indepen- 
dent states of polarization [32j . The index of refraction is 
related simply to the coherent forward-scattering ampli- 
tude / and to the number of scattering centers per unit 
volume N, in the form n = an + (2irN/k 2 ') / , where k is 
the wave number in vacuum [37j. For photons, / is the 
sum of the electronic and nuclear scattering amplitudes, 
/ = fe+fn and, for neutrons, it is the sum of the nu- 
clear and magnetic scattering lengths, / = / nuc + / mag n- 
In each homogeneous part around position r, an index of 
refraction n (r) can be defined, which can be interpreted 
as an optical potential in a wave equation 



[A + k 2 l] (r) = V (r) ip (r) 



(138) 



for a two-component wave function ip via the relation 
V (r) = 2fc 2 [(To — n (r)] [3l|. Here, the two components 
describe the two polarizations of the photon field, but 
actually a same type of equation can be used for neutrons 
as well. Indeed, based on the equation (|138|) . the elastic 
scattering of slow neutrons and of X-rays on stratified 
media have a common description (38j , which also covers 
Mossbauer scattering [3^] and diffuse scattering. The 
common scattering theory is that of the (time dependent) 
Schrodingcr equation whose stationary scattering states 
satisfy 1(155)) . 

We note that the elements of the matrix n for both 
slow neutrons and X-rays differ only slightly (typically 
10~ 5 ) from that of the 2x2 unit matrix an, which dif- 
ference can be three orders of magnitude greater, typ- 
ically 10~ 2 , in the case of Mossbauer scattering at the 
resonance energies. This difference is the basis of speak- 
ing about stronger scattering in case of the Mossbauer 
medium. In the so-defined optical medium, namely, in 
the Mossbauer medium, the magneto-optic Faraday ef- 
fect was identified [321 ) , which is often cited as a property 
of nonreciprocal media [23|. Both aspects, the presence 
of strong scattering and the magneto-optic effects, rec- 
ommend using Mossbauer medium as a model system for 
studying reciprocity. 

In the following two Sections, some examples of 
Mossbauer forward scattering on pure a— 5 Fe absorbers 
follow, to demonstrate reciprocal and nonreciprocal (in 
other words, reciprocity violating) cases. All the consid- 
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ered iron foils had the thickness of 4 /jm, and the spec- 
tra were simulated by the computer program EFFI [50jj , 
which reproduces experimental results to high precise- 
ness. As we will see, the classification of situations into 
reciprocal and nonreciprocal provides unexpected out- 
comes sometimes. 



4.1. When nonreciprocity is expected but 
reciprocity appears, instead 

It was already mentioned that, in the literature, 
magneto-optic media is often cited as nonreciprocal me- 
dia. In an early work on the field 15 1| , it was pointed out 



In an early work on the field 
that the sign of the k-parallel component of the mag- 
netization can be determined using circularly polarized 
radiation. Simply put, the situations when the (homoge- 
neous) magnetization Bjjf is parallel to the momentum 
vector k can be distinguished from the antiparallel case. 

At first sight, it seems that the here-defined arrange- 
ment is nonreciprocal. Indeed, the change of the mag- 
netization from parallel to antiparallel case can be asked 
whether it can be a 'reciprocal-plus-rotational' transfor- 
mation, discussed in Sect. [3~5l and by (|113p . the answer 
seems to be negative. In fact, a mere change of the direc- 
tion of magnetization does modify the Mossbauer spec- 
tra, as shown in Figs. [2] (a) and[2](b). 

However, taking into account the incoming and outgo- 
ing polarizations as well, if we prescribe a simultaneous 
change of the incoming and outgoing - say, right circular 
- polarizations to opposite, then the spectrum remains 
invariant, as applying (|1 13[) gives, and comparing Figs. [2] 
(a) and [H (c) justifies. Consequently, the magneto-optic 
medium discussed here itself is reciprocal (!), and this 
property either remains hidden or is revealed depending 
on how one chooses the polarizations for the processes. 

This reciprocal property can be easily understood con- 
sidering the scattering potential in the case of longitudi- 
nal Zeeman effect, in which case the potential is diagonal 
on the circular basis, fulfilling thus the requirement of 
self-transposeness, which is a manifest form of the reci- 
procity condition. 



4.2. Real nonreciprocity in Mossbauer scattering 

Unlike in the previous example, where circularly po- 
larized photons scattered on a single a — 57 Fe Mossbauer 
absorber, next considered is the case of scattering of lin- 
early polarized photons - the scattering from incident TE 
(transversal electric, in other words, er-polarization) to 
the same type of outgoing polarization - on two a — 57 Fe 
Mossbauer absorbers. The case of reciprocity for more 
then one scatterer was considered in section [331 with the 
conclusion that piecewise univectorial potentials admit a 
reciprocity operator if and only if the corresponding real 
vectors b/ of the piecewise constant potential values Vi 
are within a common plane. Since the potential in case 



of magnetic hyperfine interaction is of univectorial kind 
(see the next paragraph), in the two-scatterers case b/ 
(/ = 1,2) are necessarily within a common plane, hence, 
there are always reciprocal situations. We note that this 
is not true for the case of three scatterers with three lin- 
early independent vectors, as three such vectors do not 
lie in a common plane. 

The polarization-dependent part of the scattering am- 
plitude for Mossbauer transition in case of 57 Fe can be 
given following Eq. 5 of Ref. 52] and using the deduction 
of Rose [53j]. One finds that the scattering amplitude /, 
consequently the optical potential V, is written as a prod- 
uct of a complex number (the energy-dependent complex 
Lorentzian) and a self-adjoint matrix (the polarization- 
dependent part of /) and is, therefore, univectorial. The 
corresponding real vector b of the +1 transition is, after 
some simple algebra, found to be 



- 5 sin 2 9 sin 2ip 

-2cos6> 
sin 2 8 cos 2<p 



(139) 



where a linear basis for the polarization is used with 9 
and ip being the polar angles describing the direction of 
the hyperfine magnetic field in the laboratory system, 
which is selected so that the z axis points towards the 
propagation of the radiation and the two other axes are 
arbitrary orthogonal directions. 

Let us consider the arrangement of two samples, de- 
fined by the polar angles of the hyperfine fields 6\ — 90°, 
ipi = 90° and 6» 2 = 135°, ip 2 = 0° . The vectors b u b 2 of 
the corresponding two potentials of the two samples read 
b : = (0,0,-1) and b 2 = (0,^2, 1/2). The reciprocal 
arrangement can be defined according to the geometrical 
criterion of Sect. 13.31 namely, via the unitary transforma- 
tion that expresses a rotation around the intersection of 
the 0i— (T3 plane and the plane spanned by the vectors bi 
and b 2 , and brings the former plane to the latter. This 
arrangement has the speciality that the vector bi also 
lies in the o~\-o-$ plane and, as a consequence, it is paral- 
lel to the rotational axis, which can only be the 03 axis 
in this case. In Sect. 13.71 it was mentioned that any ro- 
tation around the 03 axis causes magnitude reciprocity. 
Indeed, Fig. [3] (a) shows that the corresponding simu- 
lated Mossbauer spectra are identical for the normal and 
reciprocal cases. 

Next, we present an arrangement, a slight modifica- 
tion of the previous example, where real observable non- 
reciprocity appears. The modification is only that we 
rotate the first foil by 45° around the axis being parallel 
to the direction of the wave propagation k, resulting in 
the polar angles of X = 90°, ip x = 45° and 6 2 = 135°, 
(f2 = 0° . Repeating the previous procedure, the b vec- 
tors of the corresponding potential read bi = (—1,0,0) 
and b 2 = (0, v2 f 1/2). In this case, bi is a vector both of 
the CT1-03 plane and of the plane spanned by the vectors 
hi and b 2 , therefore, the axis of the rotation U is the U\ 
axis. Unlike in the previous example, such a rotation does 
not cause magnitude reciprocity, thus observable nonrc- 
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ciprocity is expected. Indeed, the simulations shown in 
Fig. [3] (b) confirm the nonreciprocal property of the ar- 
rangement. Concentrating on the first Mossbauer lines 
(left ones) of the normal and reciprocal arrangements, 
one can conclude that the type of asymmetry caused by 
multiple scattering appears on the opposite side of the 
lines in the normal and the reciprocal arrangements. This 
feature fades away for thin layers (thicknesses less than 
1 /im), where the 1 st Born approximation becomes valid 
and, according to Sect. 13.91 nonreciprocity disappears. 

We note that the here-presented thought experiment 
uses the so-called polarizer-analyzer setup, therefore, if 
using radioactive source, one would lose most part of 
the intensity. A similar experiment is, however, feasible 
by using synchrotron radiation source, where the beam 
is collimated and the polarizer-analyzer setup means no 
problem. The arrangement presented in this example is 
interesting because a simple rotation by 45° switches the 
system from magnitude reciprocal to nonreciprocal. 

4.3. On the symmetry of diffuse cj-scans 

For studying lateral inhomogeneities - structural 
roughness, magnetic domains, etc. - in stratified me- 
dia, diffuse scattering, i.e., off-specular neutron [54| . 
soft- X-ray resonant magnetic diffuse scattering [55| and 
Mossbauer reflectometry (56[, is applied. A possible ex- 
perimental realization of the off-specular reflectometry is 
the so-called 'w-scan' geometry, where the detector po- 
sition is set to 29 and the sample orientation uj on the 
goniometer is varied with the sample normal remaining 
in the scattering plane (see Fig. [4]). It is straightforward 
to see in this setup that the incoming and outgoing waves 
make an angle of uj and 29 — uj with the surface of the 
stratified media, respectively. In the 'w-scan' experiment, 
the scattered intensity J as a function of u> is detected. 
In the special case uj = 9, one detects the specular ra- 
diation. A well-known property of the 'w-scan' intensity 
function is its symmetricity with respect to the specu- 
lar position, i.e., J (uj) — J (29 — uj), which is explained 
in the literature as a straightforward consequence of the 
reciprocity theorem, which is often confused with time 
reversal symmetry [571 ] . Let us now investigate how our 
reciprocity formulae can be utilized for this symmetricity 
property. 

What we ask is whether the sample at position uj, i.e., 
potential V, is a reciprocal partner of the sample at po- 
sition 29 — uj, i.e., of potential V$ — R^VR^ 1 , where 
Rrf, is the Hilbert space representation of the rotation by 
4> := (29 — uj) — uj = 2 (9 — uj) that connects the two posi- 
tions. Since the incoming and outgoing momenta are the 
same for the two arrangements, we must consider reci- 
procity combined with a rotation as done in Sect. [3~51 at 
Eq. (|116p . which we repeat here for convenience: 

(kp,pp\T\k a ,p a ) = (k p ,U R Up* a \T R \k a ,U R Up}), (140) 
where the rhs describes a scattering on RVR^ 1 = 



RUV T U~ 1 R~ 1 , with R being the Hilbert space repre- 
sentation of the rotation performing k a -R- — kp. 

The question is whether this, the scattering amplitude 
at position uj, can coincide with the one at position 29— uj, 

{kp,p 2 \T^\k a , Pl ), (141) 

with some appropriate polarizations pi, pi- To ensure 
this, on one hand we require = RVR" 1 , which can be 
expanded and rearranged as 

R^VR^ 1 = RUV T U- 1 R-\ 
R-^R^VR-^R = UV T U-\ 

R m ,,VR m l = UV T U-\ (142) 

here, in the last line, we recognized that the combination 
of the two rotations R = R^ 1 and R^ is the rotation 
R m % = Rm.K °f the sample in position uj around its 
normal m by u = 180°. Expressed in the Poincare vector 
description, condition (|142l) says 

v (O m>7r r) = v (r) (143) 
O m>7r v (O m>7r r) = Oc/Pisv (r) , (144) 

where Ou is the rotation belonging to U by We 
note that the r-dependence aspect may be successfully 
treated for lateral inhomogeneities via the DWBA ap- 
proximation. In parallel, (|144[) can be evaluated analo- 
gously to our previous analyses. 

The other requirement is that the polarizations also 
agree. Namely, if V satisfies (| 143)) — (|144|) with some U 
then one needs 

pi = U R Up*p, P 2 = U R Up* a . (145) 

Examples for such polarization settings were considered 
in Sect. E3J 

As for simple examples, the simplest one is that of a po- 
larization independent potential, V(r) — wo( r )fo- Even 
for such potentials, (| 143[) prescribes a nontrivial condi- 
tion so even polarization independent potentials must 
obey such a rotation invariance so as to exhibit the sym- 
metry of the w-scan spectrum. 

The second example is neutron scattering on a sample 
of one layer with depth independent magnetic field, which 
is a homogeneous Bj in the left half of the sample and 
a homogeneous B r in the right half (let the shape of the 
sample be left-right reflection invariant, with respect to 
a plane orthogonal to the layer). Then (|144l) requires 

Bi = O^Oc/PiaBr. (146) 

The combined transformation O m!7r O[/Pi3 is an orthog- 
onal, thus length-preserving, one. If |Bi| = |B r | then one 
can find such a definition of the o\ and 03 directions - 
the x and z directions - that (|146[) is satisfied. On the 
other hand, if |Bi| 7^ |B r | then it is impossible to fulfill 
this requirement and the spectrum of the w-scan cannot 
be symmetric. 
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5. CONCLUSION 

As demonstrated in the general-level discussion, reci- 
procity is the property of a system - describing linear 
wavelike propagation - when it is connected with the 
adjoint system by an antiunitary operator. If the reci- 
procity property is fulfilled via such a reciprocity opera- 
tor K then a reciprocity theorem holds, which expresses 
the equality of any scattering amplitude {(3\T\a) to an- 
other scattering amplitude, namely, to (Ka\T\Kf3). This 
left-right interchange of incoming and outgoing states 
gives, in important experimental applications, that a 
scattering amplitude is related to another one where the 
source and the detector are interchanged. 

Remarkably, reciprocity can hold for non-selfadjoint 
systems (systems with absorption), too. This immedi- 
ately distinguishes reciprocity from time reversal invari- 
ance. Waves with spin/polarization degree of freedom 
demonstrate that reciprocity also differs from rotational 
invariance: Rotations are unable to map an incoming po- 
larization degree of freedom to an outgoing one, nor an 
outgoing polarization to an incoming one. The above- 
presented calculations show in detail the relationship of 
reciprocity to time reversal as well as how rotation can be 
combined with reciprocity to obtain a version of the reci- 
procity theorem that is especially suitable for scattering 



experiments. 

To find reciprocity operators for a given system is a 
delicate problem, which is solved here for an important 
class of physical situations, which cover applications in 
neutron and photon scattering on multilayer structures. 
Reciprocity violation is also quantified, and the results 
are illustrated and applied on examples, chosen from the 
area of Mossbauer scattering, where reciprocity is fulfilled 
for certain processes and is immensely violated for some 
others (scattering amplitudes differing remarkably). 

The relationship established here to a recently devel- 
oping area of mathematics is expected to give an impetus 
to finding reciprocity operators to more physical systems, 
resulting in valuable applications. 
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FIG. 1: Scattering experiment in the original (upper) and 
the reversed (lower) arrangements. The reciprocity principle 
means invariance of the scattering transition amplitude under 
the reversal of the source and the detector. 
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FIG. 2: Simulated Mossbauer forward scattering spectra on 
a —° 7 Fe foil of thickness of 4 /im using the polarizer-analyzer 
setup measured for the cases: (a) incident right circular po- 
larized photons scattered to the same polarization (H > + 

scattering) and the hyperfine magnetic field Bnf being paral- 
lel to the momentum vector k, (b) incident right circular po- 
larized photons scattered to the same polarization (H > + 

scattering) and the hyperfine magnetic field But being an- 
tiparallel to the momentum vector k, and (c) incident left 
circular polarized photons scattered to the same polarization 
(— — ► — scattering) and the hyperfine magnetic field BHf 
being antiparallel to the momentum vector k. 
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FIG. 3: Simulated Mossbauer forward scattering spectra on 
two a — 57 Fc foils of thicknesses of 4 fj,m using the polarizer- 
analyzer setup for incident a (viz. transversal electric) linearly 
polarized photons scattered to the same polarization (a — > a 
scattering) for the cases the hyperfine magnetic fields Bnf,i 
and BHf,2 in the foils pointing to the directions given by the 
polar angles (a) 6i = 90°, ipi = 90° and 6 2 = 135°, <p 2 = 0°, 
and (b) 0i = 90°, ifi! = 45° and 6 2 = 135°, Lp 2 = 0°. In both 
parts, (a) and (b), solid lines denote the original and dashed 
lines the reversed (source-detector exchanged) situations. 





FIG. 4: Geometrical arrangement of the diffuse w-scan scat- 
tering experiments with 26 being the angle of the incoming 
and outgoing waves, and uj being the angle of the incident 
wave and the surface of the stratified media. The lower part 
shows the setup obtained by a rotation by <f> — 2 (9 — oS) 
around the axis being perpendicular to the scattering plane. 



